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Abstract 

The gold-plated discovery mode of a Standard Model like Higgs boson at the 
CERN Large Hadron Collider (LHC) is the H Z°Z° decay mode. To find 
and then measure the properties of the Higgs, it is crucial to have the most 
precise theoretical prediction both for the signal and the QCD background in 
this mode. In this work we calculate the effects of the initial-state multiple 
soft-gluon emission on the kinematic distributions of Z° boson and photon 
pairs produced in hadron collisions. The Collins-Soper-Sterman formalism 
is extended to resum the large logarithmic terms due to soft-gluons. The re- 
summed total rates, the invariant mass, transverse momentum, and rapidity 
distributions of the Z° and photon pairs, and the transverse momentum dis- 
tributions of the individual vector bosons are presented and compared to the 
fixed order predictions in the whole kinematic range, for the LHC energies and 
for the upgraded Fermilab Tevatron. Our conclusion is that the resummed 
predictions should be used when extracting the Higgs signal at the LHC. 
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I. INTRODUCTION 



The underlying dynamics of the electroweak symmetry breaking sector of the Standard 
Model (SM) awaits understanding. The principal goal of the CERN Large Hadron Collider 
(LHC) is to shed light on this open question. The direct searches at the CERN Large 
Electron Positron (LEP) collider have constrained the mass of the SM Higgs boson to be 
higher than 90 GeV [1]. Furthermore, global analyses of electroweak data [2] and the values 
of the top quark and the W ± boson masses [3] suggest that the SM Higgs boson is light, less 
than a few hundred GeV. Arguments based on supersymmetry (SUSY) also indicate that 
the lightest Higgs boson is lighter than the top quark [4]. Hence, the existence of a light 
Higgs boson is highly possible. 

It has been shown in the literature that a SM like Higgs boson with a mass less than 
or about 120 GeV can be detected at the upgraded Fermilab Tevatron via pp — > W ± (~^ 
£ ± u) H(-> bb,T + T-) [5], or pp(gg) -> if(-> W*W* -> £ujj and ivtv) [6], and a SUSY 
Higgs boson can be detected in the W ± h and hbb modes [7]. To observe a light (m H < 120 
GeV) neutral Higgs boson at the LHC, the most promising detection mode is the di-photon 
channel H — > 77 [8] via the production process pp(gg) — > HX. In the intermediate mass 
(120 GeV < m H < 2 m z ) region, the Z°*Z° channel is also useful in addition to the 77 
channel [8]. If the Higgs boson is heavier than twice the mass of the Z° boson, the gold- 
plated decay mode into two Z° bosons (which sequentially decay into leptons) [8] is the 
best way to detect it. At the LHC, as at any hadron-hadron collider, initial-state radiative 
corrections from Quantum Chromodynamics (QCD) interaction to electroweak processes can 
be large. Some fixed order QCD corrections have been calculated to the Higgs signal and to 
its most important backgrounds [9,10,12,13]. The next-to-leading order (NLO) corrections 
to the total cross section of pp(gg) — > HX have been found to be large (50-100 %) [9], 
and the largest contribution in the fixed order corrections results from soft gluon emission 
[14]. This signals the slow convergence of the perturbative series, and the importance of 
still higher order corrections. Furthermore, the fixed order corrections fail to predict the 
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transverse momentum distributions of the Higgs boson and its decay products correctly. 
The knowledge of these distributions is necessary to precisely predict the signal and the 
background in the presence of various kinematic cuts, in order to deduce the accurate event 
rates to compare with theory predictions [15]. To predict the correct distribution of the 
transverse momentum of the photon or Z° pair and the individual vector bosons, or the 
kinematical correlation of the two vector bosons produced at hadron colliders, it is necessary 
to include the effects from the initial-state multiple soft-gluon emission. In this work, we 
present the results of our calculation for the most important continuum backgrounds to the 
Higgs boson signal detected at hadron colliders. The distributions of the Higgs boson signal 
for the hih 2 — > H — > 77X and Z°Z°X processes, including the soft-gluon effects, will be 
discussed in our future work [16]. 

In scattering processes involving hadrons, the dynamics of the multiple soft-gluon radia- 
tion can be described by the resummation technique. We extend the Collins-Soper-Sterman 
(CSS) resummation formalism [17-19] to describe the production of photon and Z° pairs. 
This extension is analogous to our recent resummed calculation of the hadronic production 
of photon pairs [20]. In comparison, an earlier work [21] on the soft-gluon resummation 
for the qq — > Z°Z°X process did not include the complete NLO corrections. In the present 
work, the effect of initial-state multiple soft-gluon emission in qq — > Z°Z°X is resummed 
with the inclusion of the full NLO contributions, so that the inclusive rate of the Z° boson 
pair production agrees with the NLO result presented in Ref. [10]. Furthermore, we also 
include part of the higher order contributions in our results by using the CSS resummation 
formalism. 

The collected di-photon data at the Tevatron, a pp collider with center-of-mass energy 
\fS = 1.8 TeV with 84 and 81 pb -1 integrated luminosity (for CDF and DO ), are in the 
order of 10 3 events per experiment [22,23]. After the upgrade of the Tevatron, with \fS = 2 
TeV and a 2 fb" 1 integrated luminosity, about 4 x 10 4 photon pairs can be detected, and 
more than 3 x 10 3 Z° boson pairs can be produced. At the LHC, a v^S* = 14 TeV pp collider 
with a 100 fb _1 integrated luminosity, we expect about 6 x 10 6 photon and 1.5 x 10 6 Z° 
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pairs to be produced, after imposing the kinematic cuts described later in the text. This 
large data sample will play an important role in the search for the Higgs boson(s) and new 
physics that modifies the production of the vector boson pairs (e.g., by altering the vector 
boson tri-linear couplings [24]). 

The rest of this paper is organized as follows. Section II briefly summarizes the extension 
of the CSS resummation formalism to the Z°Z° pair production. In Section III, the numerical 
results of the resummed and fixed order calculations are compared for various distributions 
of the photon and Z° pairs produced at the LHC and the upgraded Tevatron. Section IV 
contains our conclusions. 



When QCD corrections to the Z° boson pair production cross section are calculated 
order by order in the strong coupling constant a s , the emission of potentially soft gluons 
spoils the convergence of the perturbative series for small transverse momenta (Qt) of the 
Z° boson pair. In the Qt <C Q region, the cross section can be written as [19] 



where Q is the invariant mass of the Z° boson pair, and the coefficients n v m are perturbatively 
calculable. At each order of the strong coupling the emitted gluon(s) can be soft and/or 
collinear, which yields a small Qt- When the two scales Q and Qt are very different, 
the logarithmic terms \n m (Q 2 /Q T ) are large, and for Q T <C Q the perturbative series is 
dominated by these terms. It was shown in Refs. [17-19] that these logarithmic contributions 
can be summed up to all order in a s , resulting in a well behaved cross section in the full Qt 
region. 

The resummed differential cross section of the Z° boson pair production in hadron colli- 
sions is written, similarly to the cross sections of the lepton pair production [25] , or photon 



II. ANALYTICAL RESULTS 



A. The CSS Resummation Formalism for Z Pair Production 
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pair production [20], in the form: 
da(h ± h 2 -> Z°Z°X) 



P 



dQ 2 dy dQ 2 T d cos 9 dcj> 48nS Q 2 



+ Y(Q T ,Q,x 1 ,x 2 ,e,<l),C 4 )}. 



(1) 



In this case, the variables Q, y, and Qt are the invariant mass, rapidity, and transverse 
momentum of the Z° boson pair in the laboratory frame, while and <fi are the polar and 
azimuthal angle of one of the Z Q bosons in the Collins-Soper frame [26]. The factor 
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originates from the phase space for producing the massive Z° boson pair. The parton 
momentum fractions are defined as X\ = e y Q/\/S, and x 2 = e~ y Q/\/~S, and \/~S is the 
center-of-mass (CM) energy of the hadrons hi and h 2 . 

The renormalization group invariant function Wij(b) sums the large logarithmic terms 
o?l \n m (b 2 Q 2 ) to all orders in a s . For a scattering process initiated by the partons i and j, 

Wij(b, Q, xi, x 2 , 9, 0, 61,2,3) = exp {-Sij(b, Q, C lj2 )} 
x \c i/hl (x 1: b, (71,2,3, t, u) C j/h2 (x 2 , b, Ci >2 , 3 , *, u) 
+ Cj/h! (xi, b, Ci )2 ,3, u, t) C i/h2 (x 2 , b, Ci j2 ,3, m, 
x^a(C 2 Q),a s (C 2 Q),0,0)- 



(2) 



Here the Sudakov exponent Q, ^1,2) is defined as 



Sij(b,Q,Ci i2 ) = J 



ciQ 2 dfi 2 



Ay (a s (/i),Ci)ln 



+ B -(a a (/Z),Ci,C 2 ) 



(3) 



In Eq. (2), JF^- originates from the hard scattering process, and will be given later for specific 
initial state partons. Ci/h{x) denotes the convolution of the perturbative Wilson coefficient 
functions C ia with parton distribution functions (PDF) f a /h(0 (describing the probability 
density of parton a inside hadron h with momentum fraction £): 



C i/h (x,b,C lt2 ,3,t,u) = 

E [ J Cia (|, b,iM = ^ , d, C 2 , t, / a/h (e, f = y ) ■ (4) 

The invariants s, t and -u are denned for the q{pi)qij>2) Z°(p 3 )Z°(p4) subprocess as 

s = (pi + p 2 ) 2 , * = (pi - Ps) 2 , w = (P2 - Ps) 2 , (5) 

with s + £ + u = 2m|. 

The functions A^, and Cjj are calculated perturbatively in powers of a s /n: 

Aj Ci) 

(« s (/i),Ci,C2) 

Cij (z,b,n,C 1 ,C 2 ,t,u) 

The dimensionless constants Ci, C 2 and C 3 = /ift were introduced in the solution of the 
renormalization group equations for Wij. Their canonical choice is C\ = C 3 = 2e~ 7B = b , 
C 2 = C1/&0 = 1, and C 4 = C 2 = 1 [19], where 7^ is the Euler constant. 

For large 6, which is relevant for small Q T , the perturbative evaluation of Eq. (2) is 
questionable. Thus in Eq. (1), is evaluated at 5* = b/^Jl + (b/b max ) 2 , so that the 
perturbative calculation is reliable. Here 6 max is a free parameter of the formalism [19] 
that has to be constrained by other data (e.g. Drell-Yan), along with the non-perturbative 
function W^ p (b) which is introduced in Eq. (1) to parametrize the incalculable long distance 
effects. Since the qq — > 77 or Z°Z°, and the qq — > V — > ££' processes have the same 
initial state as well as the same QCD color structure, in this work we assume that the 
non-perturbative function W^ p (b), fitted to existing low energy Drell-Yan data [27], also 
describes the non-perturbative effects in the qq — > 77 and Z°Z° processes. Needless to say, 
this assumption has to be tested by experimental data. 

The function Y in Eq. (1) contains contributions from the NLO calculation that are 
less singular than \jQ\ or ln(<5 2 / 'Qt) I Qt as Qt ~ * 0- This function restores the regular 
contribution in the fixed order perturbative calculation that is not included in the resummed 
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= ±(^P) n Blf(C 1 ,C 1 ), 

n=0 V 71 / 



piece Wij. In the Y function, both the factorization and the renormalization scales are chosen 
to be C4Q. The detailed description of the matching (or "switching") between the resummed 
and the fixed order cross sections for Q T ~ Q can be found in Ref. [25]. 

B. The qq, qg and qg — > Z°Z°X subprocesses 

The largest background to the Higgs boson signal in the Z°Z° channel is the continuum 
production of Z° boson pairs via the qq — > Z°Z°X partonic subprocess [28]. The next-to- 
leading order calculations of this process are given in Refs. [10,12]. A representative set of 
Feynman diagrams, included in the NLO calculations, is shown in Fig. 1. The application 
of the CSS resummation formalism for the qq — > Z°Z°X subprocess is the same as that for 
the case of qq — > 77X [20]. The and coefficients in the Sudakov exponent are 

identical to those of the Drell-Yan case. This follows from the observation that to produce 
a heavy Z° boson pair, the virtual-quark line connecting the two Z° bosons in Fig. 1 is far 
off the mass shell, and the leading logarithms due to soft gluon emission beyond the leading 
order can only be generated from the diagrams in which soft gluons are connected to the 
incoming (anti-) quark. This situation was described in more detail for di-photon production 
[20]. 

The resummed cross section is given by Eq. (1), with i and j representing quark and 
anti-quark flavors, respectively, and 

Fij(9, 9s, 0, 0) = 2Sij(gl + 9rY f~~2^ - 

The left- and right-handed couplings gL,R are defined through the qqZ° vertex, which is 
written as ry^ [g L (l - 75) + g R (l + 75)], with 

T 3 - sjQf s 2 w Q f 

9L = g and g R = -g (6) 

z,c w £c w 

Here g is the weak coupling constant, T 3 is the third component of the SU(2) L generator 
(T 3 = 1/2 for the up quark Q u , and —1/2 for the down quark Q d ), s w (c w ) is the sine (cosine) 
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of the weak mixing angle, and Qf is the electric charge of the incoming quark in the units 
of the positron charge (Q u = 2/3 and Q d = —1/3). The values of these parameters will be 
given in the next section. 

The explicit forms of the A and B functions, used in the numerical calculations are: 



A^qq (Ci) — C F , 

A^qqiCl) = C F 



Bq^iCi, C 2 ) — Cp 



(7) 



where Nf is the number of light quark flavors, N c = 3 is the number of colors in QCD, 
C F = 4/3, and ft = (11N C - 2JV / )/12. 

To obtain the value of the total cross section to NLO, it is necessary to include the Wilson 
coefficients C\f and C\f , which can be derived similarly to those for di-photon production 



[20]. The results are: 

Cjk(z,b,H, yr,t,u) = 5 jk 5(l - z), 



C 

Cjc{z,b,fi, -±,t,u) = 0, 



C { $ (z, b, ^, t, u) = 5 3k C F \Ul - z) - — In ( ^ 1 M 



1 



C${z, b, /x, t, u) = -z{\ -z)-\n[^j P^ G (z). 

In the above expressions, the splitting kernels are [29] 

'l+z 2 \ 



Ci 



+6(1 -z) 



C F 

In 2 f S^e-^ 1 - - -- -- 



D (i) 



V(t,u) , _9_ 
16 



(8) 



d(i) 

j<— A: 



» = ft 



f? M = 5 



1-3 



^ + (1 - z) 



and 



(9) 



For Z° boson pair production the function V in Eq. (7) is given by 

(-|Mf or T + 4AT(F(t, u) + F(u, t 
8 



7T 



V(t,u) = V z o z o(t,u) = — + \M 



Born | -2 




3. 



The squares of the invariant amplitudes Mq OT11 and M\} ma are defined by Eqs.(6) and (7) of 
Ref. [11]. For Z° pair production N = N c Q A e+ ^ e ((g1 Zq ) 4 + (g q + q ) 4 ), where N c = 3 is the 
number of colors in QCD, Q e + is the electric charge of the positron, 4 — 2e is the number 
of space-time dimensions, and g*± q are defined by Eq.(3) of Ref. [11]. The definition of 
the function F(t,u) is somewhat lengthy and can be found in Appendix C of Ref. [10] (cf. 
Eqs. (CI) and (C2)). The function V(t,u) depends on the kinematic correlation between 
the initial and final states through its dependence on t and u. In the mz — > limit, 
F(t,u) + F(u,t) reduces to F mrt (t,u) of the di-photon case which is given in Ref. [20]. 1 
The non-perturbative function used in this study is [27] 

w Si P (p } Q, Qo, x 1: x 2 ) = exp 

with 9l = 0.11 GeV 2 , g 2 = 0.58 GeV 2 , g 3 = -1.5 GeV _1 , and Q = 1.6 GeV. These values 
were fit for the CTEQ2M parton distribution function, with the canonical choice of the 
renormalization constants, i.e. C\ = C3 = 60 and C2 = 1, and b max = 0.5 GeV -1 was used. 
In principle, these coefficients should be refit for CTEQ4M distributions [31] used in this 
study. We have checked that using the updated fit in Ref. [32] does not change largely our 
conclusion because at the LHC and Tevatron energies the perturbative Sudakov contribution 
is more important compared to that in the low energy fixed target experiments. 

Before concluding this section we note that for the di-photon production, we use the 
formalism described in Ref. [20] to include the gg — > 77X contribution, in which part of 
the higher order corrections has been included via resummation. Since a gauge invariant 
calculation of the gg — > Z°Z°g cross section in the SM involves diagrams with the Higgs 
particle, we shall defer its discussion to a separate work [16]. 

1 This is connected to the fact that as mz — ► the virtual corrections of the Z° pair and di- 
photon productions are the same (up to the couplings), which is apparent when comparing Eq.(ll) 
of Ref. [30] and Eq.(12) of Ref. [10], after including a missing factor of l/(167rs) in the latter 
equation. 
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III. NUMERICAL RESULTS 



We implemented our analytic results in the ResBos Monte Carlo event generator [25]. 
As an input we use the following electroweak parameters [33]: 

Gp = 1.16639 x 1CT 5 GeV~ 2 , m z = 91.187 GeV, m w = 80.41 GeV, a(m z ) 1 



128.88 

In the on-shell renormalization scheme we define the effective weak mixing angle 



,2 



with 



m 



P = 



w 



Tflr 



w 2 ' 



f na{m z ) \ 1 
v V2G F m 2 w ) 



'■z 

In Eq.(6), the coupling of the Z° boson to fermions, g, is defined using the improved 
Born approximation: 

g 2 = 4V2G F (c^) 2 m 2 zP , 

with c e JJ — \ 1 — sin e e J f , the cosine of the effective weak mixing angle. (In Eq. (6) c w is 
identified with c e J f .) We use the NLO expression for the running strong and electroweak cou- 
plings a s (fj) and a(/i), as well as the NLO parton distribution function CTEQ4M (defined in 
the modified minimal subtraction, i.e. MS, scheme), unless stated otherwise. Furthermore, 
in all cases we set the renormalization scale equal to the factorization scale: /ir — /if = Q- 
Table I summarizes the total rates for the leading order (LO), i.e. 0(a®), and resummed 
photon- and Z°-pair production cross sections for the LHC and the Tevatron. For the 
lowest order calculation we show results using LO (CTEQ4L) and NLO (CTEQ4M) parton 
distributions, because there is a noticeable difference due to the PDF choice. As it was 
discussed in Ref. [25], the resummed total rate is expected to reproduce the 0(a s ) rate, 
provided that in the resummed calculation the and coefficients and the 0(a s ) 

Y piece are included, and the Qt distribution is described by the resummed result for 
Qt < Q and by the 0(a s ) result for Q T > Q. In our present calculation we added the 
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coefficient to include the most important higher order corrections in the Sudakov exponent. 
Our matching prescription (cf. Ref. [25]) is to switch from the resummed prediction to the 
fixed-order perturbative calculation as they cross around Qt ~ Q. This switch is performed 
for any given Q and y of the photon or Z° boson pairs. In the end, the total cross section 
predicted by our resummed calculation is about the same as that predicted by the NLO 
calculation. The small difference of those two predictions can be interpreted as an estimate 
of the contribution beyond the NLO. 

A. Z° pair production at the LHC 

In the LHC experiments the H — > zJ°zJ° channel can be identified through the decay 
products of the Z° bosons. The detailed experimental kinematic cuts for this process are 
given in Ref. [8]. Since the aim of this work is not to analyze the decay kinematics of the 
background, rather to present the effects of the initial-state soft-gluon radiation, following 
Ref. [10] for the LHC energies we restrict the rapidities of each Z° bosons as: \y z \ < 3.0. 
We do not apply any other kinematic cuts. The total rates are given in Table I. Our 
O{a s ) rates are in agreement with that of Ref. [10] when calculated using the same PDF. 
We expect the resummed rate to be higher than the 0(a s ) rate due to the inclusion of 
the term. Indeed, our K factor, defined as the ratio of the resummed to the LO rate 
using the same PDF in both calculations, is higher than the naive soft gluon ^-factor 
(-^naive = 1 + 87ra s (Q)/9 ~ 1.3) of Ref. [34], which estimates the NLO corrections to the 
production rate of qq — > Z°Z°X in the DIS (deep-inelastic scattering) scheme. Our if-factor 
approaches the naive one with the increase of the center-of-mass energy, as expected. 

The rates for the different subprocesses of the Z° boson pair production are given in 
Table II. At the LHC the qg — > Z°Z°X subprocess contributes about 25% of the qq + qg — > 
Z°Z°X rate. The A"-factor is defined as the ratio cr(qq + qg — > Z°Z°X)/a(qq — > Z°Z°), 
which is about 1.4 for using CTEQ4M PDF. 

Figs. 2-4 show our results for proton-proton collisions at the LHC energy, v^S* = 14 GeV. 
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Fig. 2 shows the transverse momentum distribution of Z° pairs. The NLO (0(a s )) prediction 
for the qq + qg — > Z°Z°X subprocesses, shown by the dotted curve, is singular as Qt — > 0. 
This singular behavior originates from the contribution of terms which grow at least as fast 
as 1/Qt or ln(Q 2 / Q T ) / Qt ■ This, so-called asymptotic part, is shown by the dash-dotted 
curve, which coincides with the 0(a s ) distribution as Qt — > 0. After exponentiating these 
terms, the distribution is well behaved in the low Qt region, as shown by the solid curve. 
The resummed curve matches the 0(a s ) curve at about Qt = 320 GeV. Following our 
matching prescription described in the previous section, we find that this matching takes 
place around Qt = 300 GeV, depending on the actual values of Q and y. Fig. 2 also shows 
that at the LHC there is a substantial contribution from qg scattering, which is evident from 
the difference between the solid and dashed curves, where the dashed curve is the resummed 
contribution from the qq — > Z°Z°X subprocesses. 

In Fig. 3 we give the integrated distributions, defined as 



where Q™ ax is the largest Q T allowed by the phase space. In the NLO calculation, this 
distribution grows without bound near Q^ m = 0, as a result of the singular behavior of the 
scattering amplitude when Qt — > 0. It is clearly shown by Fig. 3 that the Qt distribution 
of the resummed calculation is different from that of the NLO calculation. The different 
shapes of the two curves in Fig. 3 indicates that the predicted Z° pair production rates, 
with a minimal value of the transverse momentum Qt, are different in the two calculations. 
This is important at the determination of the background for the detection of a Higgs boson 
even with moderately large transverse momentum. For Q™ in = 50 GeV, the resummed cross 
section is about 1.5 times of the NLO cross section. 

The invariant mass and the rapidity distributions of the Z° boson pairs, and the trans- 
verse momentum distribution of the individual Z° bosons are shown in Fig 4. When calcu- 
lating the production rate as the function of the Z° pair invariant mass, we integrate the 
Qt distribution for any Q, and y. When plotting the transverse momentum distributions of 




(10) 
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the individual Z° bosons, we include both of the Z° bosons per event. In the shape of the 
invariant mass and rapidity distributions we do not expect large deviations from the NLO 
results. Indeed, the shape of our invariant mass distribution agrees with that in Ref. [10]. 
However, the resummed transverse momentum distribution P z of the individual Z° bosons 
is slightly broader than the NLO distribution (not shown in Fig 4, cf. Ref. [10]). This is 
expected because, in contrast with the NLO distribution, the resummed transverse momen- 
tum distribution of the Z° boson pair is finite as Qt — > so that the P z distribution is less 
peaked. 

In Figs. 5 and 6 we show the dependence of the resummed result on the values of the 
renormalization constants C, (i = 1,2), and the values of the non-perturbative parameters 
Qi (i = 1,2,3), respectively. As Eq. (2) shows, both the Sudakov exponent and the 
Wilson coefficients Ci/h depend on the renormalization constants C\ and C 2 . The scale C\jb 
determines the onset of the non-perturbative physics, and C 2 Q specifies the scale of the 
hard scattering process. We vary both C\ and C 2 by a factor of 2. In Fig. 5, we show that 
the resummed calculation using the canonical C\ and C 2 values (upper solid curve) hardly 
differs from the one which uses C\ = b /2 and C 2 = 1/2 (lower solid curve). This difference is 
certainly smaller than the difference between the resummed and fixed order (dashed) curves 
in the Q zz = 50 — 100 GeV region. In Fig. 6, we show two resummed curves, one with the 
non-perturbative parameters given at the end of Section II. B (solid curve), and one with the 
following non-perturbative parameters: g 1 = 0.15 GeV 2 , g 2 = 0.48 GeV 2 , g 3 = —.58 GeV -1 , 
and Qo — 1.6 GeV (dashed curve, cf. Ref [32]). There is some difference only in the lowest 
Q zz region (Q zz < 30 GeV), and this difference is negligible compared to the difference 
between the resummed and NLO (dotted) calculations. Based on these results, we conclude 
that the CSS resummation gives a stable prediction for the gauge boson pair production at 
the LHC energy, and the same conclusions also hold for the Tevatron. 
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B. Z° Pair Production at the upgraded Tevatron 



After the upgrade of the Fermilab Tevatron, there are more than 3 x 10 3 Z° boson pairs to 
be produced. Since this data sample can be used to test the tri-linear gauge boson couplings 
[24], we also give our results for the upgraded Tevatron with proton-anti-proton collisions 
at a center-of-mass energy of 2 TeV. Our kinematic cuts constrain the rapidity of both of 
the Z° bosons such that \y z \ < 3. Both the LO and resummed total rates are listed in 
Table I. The ratio a(qq + qg — > Z°Z°X)/a(qq Z°Z°) is about 1.6, which is larger than 
the naive soft gluon ^-factor of 1.3. Table II shows that qg — > Z°Z°X partonic subprocess 
contributes only a small amount (about 3%) at this energy, in contrast to 25% at the LHC. 

Figs. 7-9 show the resummed predictions for the upgraded Tevatron. The invariant mass 
and rapidity distributions of Z° boson pairs, and the transverse momentum distribution 
of the individual Z° bosons are shown in Fig. 7. The solid curve shows the resummed 
contributions from the qq + qg — > Z°Z°X subprocess. The resummed contribution from 
the qq — > Z°Z°X subprocess is shown by the dashed curve. The leading order qq — > Z°Z° 
cross section is also shown, by the dash-dotted curve. The invariant mass distribution of 
the qq + qg — > Z°Z°X subprocess is in agreement with the NLO result of Ref. [10], when 
calculated for v^S* =1.8 TeV. From this figure we also find that the contribution from the 
qg — > Z°Z°X subprocess at the energy of the Tevatron is very small. 

In Fig. 8 we compare the NLO and resummed distributions of the transverse momentum 
of the Z° pair. The figure is qualitatively similar to that at the LHC, as shown in Fig. 2. 
The resummed and the NLO curves merge at about 100 GeV. The resummed contribution 
from the qq — > Z°Z°X subprocess is shown by the dashed curve, which clearly dominates 
the total rate. 

In Fig. 9 we show the integrated distribution da/dQ'™ 11 for Z° boson pair production at 
the upgraded Tevatron. The figure is qualitatively the same as that for the LHC (cf. Fig. 3). 
The NLO curve runs well under the resummed one in the Q™ m < 80 GeV region, and the 
Qt distributions from the NLO and the resummed calculations have different shapes even 
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in the region where Qt is of the order 60 GeV. For Q™ m = 30 GeV, the resummed rate is 
about 1.5 times of the NLO rate. 

C. Di-photon production at the LHC 

Photon pairs from the decay process H — > 77 can be directly detected at the LHC. When 
calculating its most important background rates, we impose the kinematic cuts on the final 
state photons that reflect the optimal detection capabilities of the ATLAS detector [8]: 

Pj. > 25 GeV, for the transverse momentum of each photons, 

|y 7 | < 2.5, for the rapidity of each photons, and 

Prl (.Pt+Pt) < 0-7) to suppress the fragmentation contribution, where p\ is the transverse 
momentum of the photon with the higher p T value. 

We also apply a AR = 0.4 separation cut on the photons, but our results are not sensitive 
to this cut. (This conclusion is similar to that in Ref. [20].) The total rates and cross 
sections from the different partonic subprocesses are presented in Tables I and III. We 
have incorporated part of the higher order contributions to this process by including in 
the Sudakov factor and in the Wilson coefficient functions (cf. Ref. [20]). Within this 
ansatz, up to O(a^), the gg — > 77X rate is about 24 pb, which increases the total if-factor 
by almost 1.0. The leading order gg — > 77 rate, via the box diagram, is about 22 pb and 14 
pb for using the LO PDF CTEQ4L and the NLO PDF CTEQ4M, respectively. The large 
difference mainly due to the differences in the strong coupling constants used in the two 
calculations 2 : CTEQ4L requires a s {m z ) = 0.132, while for CTEQ4M a s {m z ) = 0.116. The 
ratio a(qq + qg — > r yyX)/a(qq — > 77) is 1.5, and a(gg — > r yyX)/a(qq — > 77) is about 1. 
Hence, the ratio of the resummed and the O(a s ) rates, is quite substantial. 

Figs. 10-12 show our predictions for distributions of di-photons produced at the LHC. 
In Fig. 10 we plot the invariant mass and rapidity distribution of the photon pairs, and the 



2 When using the CTEQ4L PDF, we consistently use the LO running coupling constant a, 
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transverse momentum distribution of the individual photons. When plotting the transverse 
momentum distributions of the individual photons we include both photons per event. The 
total (upper solid) and the resummed qq + qg — > 77X (dashed), qq — > 77X (dotted), 
99 — > 11 X (dash-dotted), and the fragmentation (lower solid), as well as the leading order 
qq — > 77 (middle solid) contributions are shown separately. The ratio of the resummed and 
the LO distributions is about 2.5 which is consistent with the result in Table I. The relative 
values of the contributions from each subprocesses reflect the summary given in Table III. 

Fig. 11 shows various contributions to the transverse momentum of the photon pair. 
At low Qt values (Qt <C Q), the qq — > 77X contribution is larger than the qg — > 77X 
contribution, while at high Q T values (Qt > Q), the qg — * 77X subprocess becomes more 
important. The gg contribution dominates the total rate in low Qt region, and the kink in 
the gg curve at about 50 GeV indicates the need for the inclusion of the complete 0{a\) 
gg — > 77(7 contribution. (Recall that our prediction for the gg contribution at 0(a^) only 
holds for small Q T) where the effect of the initial-state soft-gluon radiation is relatively more 
important for a fixed Q.) 

In Fig. 12 we give the integrated cross section as the function of the transverse momentum 
of the photon pair produced at the LHC. Similarly to the Z° pair production, there is a 
significant shape difference between the resummed and the NLO curves in the low to mid 
Qt region. For Q^ m = 50 GeV, the resummed rate is about 1.5 times of the NLO rate. 

D. Di-photon production at the upgraded Tevatron 

In Ref. [20], we have presented the predictions of the CSS resummation formalism for 
the di-photon production at the Tevatron with v^S* = 1.8 TeV, and compared with the data 
[22,23]. In this paper, we show the results for the upgraded Tevatron with \fS = 2.0 TeV. 
We use the same kinematic cuts which were used in Ref. [20]: 

p T > 12 GeV, for the transverse momentum of each photons, 

|y 7 | < 0.9, for the rapidity of each photons. 
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An isolation cut of AR = 0.7 is also applied. The total cross sections and the rates of the 
different subprocesses are given by Tables I and III. The ratio of the qq + qg — > 77X 
and qq — > 77 rates is about 1.5, similar to that at the LHC. The leading order rate for the 
gg -> 77 subprocess is about 6.0 pb and 4.3 pb for using CTEQ4L and CTEQ4M PDF, 
respectively. The NLO rate for gg — ■> 77$ is estimated to be 8.3 pb, using the approximation 
described in Ref. [20], which is about the same magnitude as the leading order qq — > 77 rate. 
From our estimate of the NLO gg rate, we expect that the complete O(a^) contribution will 
be important for photon pair production at the Tevatron. 

Figs. 13-15 show our results for photon pairs produced at the upgraded Tevatron. The 
resummed predictions for the invariant mass and rapidity distributions of the photon pairs, 
and the transverse momentum distribution of the individual photons are shown in Fig. 13. 
In Fig. 14 we also plot the contributions to the transverse momentum of the photon pair 
from the qq + qg — > 77X (dashed), qq — > 77X (dotted), gg — > 77(7 (dash-dotted), and 
the fragmentation (lower solid) subprocesses, separately. The leading order qq — > 77 cross 
section (middle solid) is also plotted. In the low Q T region, the gg and the qq rates are about 
the same, and the qg rate becomes more important in the large Qt region. Furthermore, 
after imposing the above kinematic cuts, the fragmentation contribution is found to be 
unimportant. 

Fig. 15 shows the integrated Q T distribution. The qualitative features of these distribu- 
tions are the same as those predicted for the LHC. For Q!f in = 10 GeV, the resummed cross 
section is about twice of the NLO cross section. 

IV. CONCLUSIONS 

In this work we studied the effects of the initial-state multiple soft-gluon emission on the 
total rates and various distributions of the most important background processes (pp, pp — > 
77X, Z°Z°X) to the detection of the Higgs boson at the LHC. We applied the extended 
CSS formalism to resum the large logarithms induced by the soft-gluon radiation. We 
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found that for the qq and qg initiated processes, the total cross sections and the invariant 
mass distributions of the photon and Z° boson pairs are in agreement with the fixed order 
calculations. From our estimate of the NLO rate of the gg initiated process, we expect 
that the complete O(a^) contribution will be important for photon pair production at the 
Tevatron. We showed that the resummed and the NLO transverse momentum distributions 
of the Z° and photon pairs are substantially different for Qt ~ Q/2. In terms of the 
integrated cross section above a given <3™ m , this difference can be as large as 50% in the 
low to mid-range of Q™ m . Using the resummation calculation, we are able to give a reliable 
prediction of the Qt and any other distribution in the full kinematical region at the LHC 
and the Tevatron, even in the presence of kinematic cuts. Since the bulk of the signal is in 
the low transverse momentum region, we conclude that the difference between the NLO and 
resummed predictions of the background rates will be essential when extracting the signal 
of the Higgs boson at hadron colliders. 
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Di-boson 


Collision 


Vs 


Fixed Order O(a s ) 


Resummed 


produced 


type 


(TeV) 


CTEQ4L 


CTEQ4M 


© 0(a s ) 


Z°Z° 


pp 


14 


9.14 


10.3 


14.8 


Z°Z° 


pp 


2 


0.91 


1.01 


1.64 


77 


pp 


14 


22.1 


24.5 


60.8 


77 


pp 


2 


8.48 


9.62 


22.8 


77 


pp 


1.8 


6.30 


7.15 


17.0 



TABLE I. Total cross sections of di-photon and Z° boson pair production at the LHC and the 
upgraded Tevatron, in units of pb. The kinematic cuts are described in the text. The "©" sign 
refers to the matching prescription discussed in the text. 





Collision 


qq -► Z°Z°X 


qg -»• Z°Z°X 


(TeV) 


type 






14 


pp 


10.9 


3.91 


2 


pp 


1.62 


0.02 



TABLE II. Resummed cross sections of the subprocesses for Z° boson pair production at the 
LHC and the upgraded Tevatron, in units of pb. The kinematic cuts are described in the text. 
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Vs 


Collision 


qq -> 77X 


95 77^ 


ffff 77 


99 775 


g# -> 7gX -> 77 A" 


(TeV) 


type 






0(a 2 s ) 4L O(a^) 4M 




Fragmentation 


14 


pp 


20.5 


16.6 


22.3 14.4 


23.9 


6.76 


2 


pp 


9.68 


4.81 


6.02 4.34 


8.26 


2.15 



TABLE III. Cross sections of the subprocesses for di-photon production at the LHC and the 
upgraded Tevatron, in units of pb. The resummed qg — > 77A rate includes the fragmentation 
contribution. The 0(a 2 s ) gg — > 77 rates were calculated using both the CTEQ4L and CTEQ4M 
PDF's. The kinematic cuts are described in the text. 
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FIG. 1. A representative set of Feynman diagrams included in the NLO calculation of Z° pair 
production. 
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FIG. 2. Transverse momentum distribution of Z° pairs from qq + partonic initial states at 
the LHC. The 0(a s ) (dotted) and the asymptotic (dash-dotted) pieces are coincide and diverge 
as Qt — ► 0. The resummed (solid) curve matches the 0(a s ) curve at about Qt = 320 GeV. The 
resummed qq contribution (excluding the qg contribution) is shown as dashed line. 
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FIG. 3. The integrated cross section for Z° boson pair production at the LHC. The resummed 
and the 0(a s ) distributions are shown in solid and dashed lines, respectively. 
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FIG. 4. Invariant mass and rapidity distributions of Z° boson pairs, and transverse momen- 
tum distributions of the individual Z° bosons at the LHC. The resummed contribution of the 
QQ + 19 ~^ Z°Z°X subprocess is shown by the solid curve, and of the qq — ► Z°Z°X subprocess by 
the dashed curve. The leading order distribution of qq — > Z°Z° is shown by the dash-dotted curve. 
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FIG. 5. Resummed and NLO transverse momentum distributions of Z° boson pairs at 
LHC. The two resummed curves are calculated for C\ = bo and C2 = 1 (upper solid), and 
C\ = 6o/2 and Ci = 1/2 (lower solid), respectively. The NLO curves are the same as in Fig. 2. 
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FIG. 6. Resummed and NLO transverse momentum distributions of Z° boson pairs at the LHC. 
The two resummed curves are calculated for gi = 0.11 GeV 2 , g 2 = 0.58 GeV 2 , and g 3 = -1.5 GeV" 1 
(solid), and for and for g x = 0.15 GeV 2 , g 2 = 0.48 GeV 2 , and g 3 = -.58 GeV" 1 (dashed), 
respectively. In both cases Qo = 1.6 GeV was used. The NLO (dotted) curve is the same as in 
Fig. 2. 
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FIG. 7. Same as Fig. 4 except for the upgraded Tevatron. 
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FIG. 8. Same as Fig. 2 except for the upgraded Tevatron. 
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FIG. 9. Same as Fig. 3 except for the upgraded Tevatron. 
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FIG. 10. Invariant mass and rapidity distributions of photon pairs, and transverse momentum 
distributions of the individual photons at the LHC. The total resummed contribution (upper solid), 
and the resummed qq + qg — > 77X (dashed), qq — > 77X (dotted), gg — > 775 (dash-dotted), as 
well as the fragmentation (lower solid) contributions are shown separately. The qq — ► 77 <9(a£) 
distribution is shown in the middle solid curve. 
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FIG. 11. Transverse momentum distribution of photon pairs at the LHC. The total resummed 
contribution (upper solid), the resummed qq+qg — ► 77X (dashed), qq — > 77X (dotted), 55 — > 77X 
(dash-dotted), as well as the fragmentation (lower solid) contributions are shown separately. 
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FIG. 12. The integrated cross section for photon pair production at the LHC. The resummed 
and the 0(a s ) distributions are shown in solid and dashed lines, respectively. 
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FIG. 13. Same as Fig. 10 but for the upgraded Tevatron. 
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FIG. 14. Same as Fig. 11 but for the upgraded Tevatron. 
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FIG. 15. Same as Fig. 12 but for the upgraded Tevatron. 
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